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PURPOSE:  The  Coastal  and  Hydraulics  Engineering  Technical  Note  (CHETN)  herein 
describes  a  method  for  estimating  the  percentage  of  sand-sized  material  placed  in  suspension  by 
breaking  waves  and  carried  by  a  cross-channel  current  to  either  fall  into  a  navigation  channel  or 
travel  across  it.  Required  inputs  are  channel  width  and  depth,  depth  in  the  vicinity  of  the 
channel,  depth-averaged  current  velocity  perpendicular  to  the  channel,  and  sediment  fall  speed. 
The  procedure  is  applicable  to  inlet  entrances  that  experience  breaking  waves.  The  percentage  of 
material  deposited  into  or  passing  a  channel  is  an  input  to  the  channel  infilling  model  described 
in  Kraus  and  Larson  (2001). 

CALCULATION  METHOD:  Navigation  channels  issuing  from  an  inlet  intercept  sediment  that 
is  moving  alongshore  (Figure  1).  The  longshore  transport  may  be  generated  by  wave-  and  wind¬ 
generated  currents,  and  by  the  longshore  component  of  the  flood-tidal  current  entering  the 
channel.  Sediment  can  pass  over  the  channel  by  moving  in  suspension,  and  it  can  be  deposited  on 
the  channel  bottom  and,  possibly,  resuspended  and  transported  out  of  the  channel  (Kraus  and 
Larson  2001).  This  Technical  Note  presents  a  method  for  estimating  the  percentage  of 
suspended  material  transported  alongshore  (perpendicular  to  a  channel)  that  will  be  deposited  in 
a  channel  (the  trapping  ratio),  from  which  the  amount  passing  over  the  channel  is  determined. 
The  following  assumptions  are  made: 

a.  Transport  occurs  mainly  in  the  surf  zone,  for  example,  in  Sections  2,  3,  and  4  in  Figure  1. 
This  assumption  will  be  removed  in  later  development  of  the  model  to  describe  deep- 
draft  channels  and  transport  seaward  of  the  surf  zone. 

b.  Deposition  in  the  channel  is  controlled  by  gravity  through  the  sediment  fall  speed  Vf. 
Upward  diffusion  by  turbulence  and  vertical  currents  in  the  channel  is  neglected;  giving  a 
conservative  estimate  of  infilling. 

c.  The  flow  across  the  channel  can  be  described  by  a  simple  parameterization  based  upon 
continuity  (recirculation  in  the  channel  is  neglected). 

d.  The  velocity  is  uniform  through  the  water  column. 
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Figure  1 .  Planview  sketch  for  channel  transport  sections 


Transport  at  Channel:  The  transport  rate  qin  at  the  updrift  edge  of  the  channel  may  be 
calculated  from  the  depth-averaged  longshore  current  (current  uniform  through  the  water 
column)  and  the  suspended  sediment  concentration  distribution  according  to  the  product, 

ha  ha 

qin  =  \UaC{z)dz  =  Ua\C{z)dz  (1) 

0  0 

where  Ua  =  depth-averaged  velocity  normal  to  the  channel,  C(z)  =  sediment  concentration,  z  = 
elevation  from  the  surrounding  (ambient)  bed,  and  ha  is  the  ambient  depth  adjacent  to  the 
channel  (Figure  2).  The  transport  rate  per  unit  length  of  channel  may  have  to  be  divided  into 
sections  as  shown  in  Figure  1,  according  to  the  current,  ambient  depth,  channel  depth,  and 
presence  of  wave  breaking. 

Empirically,  the  sediment  concentration  profile  in  the  surf  zone  has  been  found  to  decrease  with 
elevation  from  the  bottom  as  an  exponential  shape, 

C(z)  =  Cb  exp (-kz !ha)  (2) 

where  C&  =  concentration  at  the  bed,  and  A  =  dimensionless  empirical  suspended-sediment  decay 
coefficient  (Kraus  and  Dean  1987).  The  transport  rate  of  sediment  suspended  to  a  certain  level  zs 
above  the  ambient  bed  adjacent  to  the  channel  is  then  found  by  integration  to  that  elevation  as: 

0, )  =  u a  f  cb  exP  ( -'kz  /  K ) dz  =  U“Cbha  [  1  -  cxp(-Xzs  !ha)\  (3) 

0  K 
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If  zs  -  ha,  the  total  transport  through  the  water  column  is  obtained: 
?»A)  =  ^YiL[l-exp(-X)] 


(4) 


The  advective  velocity  across  the  channel  Uc  can  be  estimated  from  the  continuity  equation 
given  the  velocity  and  ambient  depth  upstream  of  the  channel  and  the  average  depth  hc  in  the 
channel: 


h 

U  =-2-U 

c  h  a 


(5) 


Therefore,  if  a  suspended  sediment  particle  of  fall  speed  fy  is  transported  with  advective  velocity 
Uc  across  a  channel  of  width  W,  it  will  fall  a  vertical  distance  A z  given  by: 


W  h  Vf 
Az-  — Vf=  ——W 


U. 


h,  U„ 


(6) 


Particles  suspended  to  level  A z  or  greater  at  the  updrift  end  of  the  channel  will  pass  over  the 
channel,  whereas  particles  with  elevations  less  than  A z  will  fall  into  the  channel  and  be  deposited 
there. 


The  ratio  p  of  the  total  suspended  transport  qin  that  is  trapped  by  the  channel  can  be  estimated 
from: 
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^„(Az)_1-exP(-^Az//^) 

VJK)  I  —  cxp(— A) 


1-exp 


f  h  W  V  ^ 

-Jt-5 - L 


V  K  U.J 


1  -  exp(-A,) 


(7) 


With  the  above  expressions,  p,  which  can  be  called  a  “trapping  ratio,”  can  be  estimated  from  the 
idealized  geometry  of  the  channel  (channel  width  and  average  depth  at  the  particular  time  of  the 
calculation),  sediment  fall  speed,  upstream  current  speed  and  depth,  and  decay  coefficient  X  of 
the  concentration  profile. 

Table  1  lists  sediment  fall  speeds  for  quartz  grains  of  typical  diameters  and  for  different  water 
temperatures,  which  will  change  water  viscosity  and,  therefore,  sediment  fall  speed. 


Table  1 

Selected  Values  of  Fall  Speeds  (m/s)  for  Quartz  Sand  Particles  of  Specified  Diameters 
(after  Hallermeier  (1981)) 

Temp, 
deg  C 

Grain  Size,  mm 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.50 

0.60 

0.80 

1.00 

2.00 

10 

0.015 

0.023 

0.029 

0.035 

0.042 

0.048 

0.062 

0.075 

0.104 

0.132 

0.189 

20 

0.018 

0.025 

0.032 

0.039 

0.046 

0.054 

0.069 

0.084 

0.115 

0.133 

0.189 

30 

0.020 

0.02/ 

0.035 

0.043 

0.051 

0.059 

0.075 

0.092 

0/119 

0/133 

0.189 

Suspended  Sediment  Decay  Coefficient:  The  decay  coefficient  X  has  been  measured  to 
have  a  representative  value  of  1.65  with  standard  deviation  of  0.68  for  surf-zone  sand  with 
median  grain  size  in  the  range  of  0.14  to  0.22  mm  (Kraus  and  Dean  1987).  Because  X  is  grain- 
size  dependent,  in  this  section  a  rational  method  is  presented  for  extending  the  result  to  other 
grain  sizes. 

If  the  time-averaged  turbulence  intensity  is  homogeneous  through  the  water  column,  which  is  a 
reasonable  representation  of  a  surf  zone  or  area  of  breaking  waves,  it  is  known  that  the  vertical 
distribution  of  the  sediment  concentration  is  proportional  to  the  quantity  exp  [-(f)  /  s  s)z] ,  where 

s.5=  sediment-mixing  coefficient.  Adopting  concepts  of  Battjes  (1975),  the  sediment-mixing 
coefficient  can  be  estimated  as, 


f  £)  V/3 

Py 


(8) 


where  kci  =  dimensionless  empirical  coefficient  expected  to  be  independent  of  grain  size,  D  = 
energy  dissipation  produced  by  the  breaking  waves,  and  p  =  density  of  water.  The  standard 
expression  expf-^- /ss)z]  becomes, 


C(z)  =  Cb  exp 


Vf 


V  K(D/p)mhaJ 


(9) 
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and  comparison  of  Equation  9  with  the  empirical  Equation  2  gives  the  decay  coefficient  X  as: 


X  = 


Vf 

kd{D!  p)1/3 


(10) 


The  dissipation  D  can  be  estimated  if  the  nearshore  profile  updrift  the  channel  has  an 
approximate  equilibrium  shape.  The  development  of  Dean  (1977)  then  gives  the  energy 
dissipation  as: 


D 

P 


3/2  2  ,3/2, 

=  —g  Y  bA  K 


(11) 


where  yb  =  breaker  index  (0.8),  g  =  acceleration  due  to  gravity,  and  A  =  equilibrium  profde  shape 
parameter,  which  Kriebel,  Kraus,  and  Larson  (1991)  showed  could  be  estimated  by: 


A  = 


9 

4 


4F2V/3 
_L 

g 


(12) 


Substituting  the  expressions  for  D/p  and  A  into  X  yields: 


X  =  - 


1/3 


4L 


§K 


(13) 


Equation  (13)  determines  X  as  a  function  of  the  sediment  fall  speed  and  ambient  depth,  and  it  is 
expected  to  have  the  correct  functional  dependencies.  Quantitative  accuracy  must  be  determined 
empirically  through  the  value  of  kd,  which  is  expected  to  be  independent  of  environmental 
factors.  Based  on  the  sediment  streamer  trap  data  of  Kraus  and  Dean  (1987)  for  nominal  0.2-mm 

particle  size  sand,  kd  =  0.03.  The  quantity  Vf  /  yjgha  appearing  in  Equation  13  is  a  Froude 

number.  For  typical  sand  sizes  with  fall  speeds  in  the  range  of  0.02  to  0.06  m/s  (Table  1),  and 
typical  ambient  depths  of  1  to  10  m,  this  quantity  is  in  the  range  of  10~2  to  10"3. 

Figure  3  plots  values  of  p  (in  percent)  as  a  function  of  the  argument  hcWVf  /(h2aUa )  in  Equation  7 
for  values  of  the  Froude  number  Fr  =  V f  /  ^Jgh,,  ranging  from  0  to  0.02  in  four  steps  of  0.005 

and  kd  =0.03.  For  sand-sized  and  finer  sediments,  Fr  is  expected  to  be  small,  and  values  of p  will 
lie  close  to  a  straight  line.  A  small  value  on  Fr  implies  a  small  value  of  X  (Equation  13), 
describing  a  concentration  profile  that  is  approximately  uniform  through  the  water  column.  If  Fr 
approaches  zero  (i.e.,  k— >0),  the  trapping  ratio  p  will  be  close  to  hcWVf  /(h2aUa) ,  which 

corresponds  to  the  portion  of  the  water  column  upstream  of  the  channel  that  Az  constitutes. 
However,  because  the  quantity  hcWVf  l{h2aUa)  contains  Fr,  the  limit  Fr  =  0  implies  that  this 

quantity  is  zero  as  well.  Thus,  in  Figure  3,  the  line  labeled  “0+”  denotes  an  Fr- value  slightly 
greater  than  zero.  In  this  limit,  for  which  Fr  =  0+  (i.e.,  A,=0),  the  trapping  ratio  will  be 
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p  =  hcWVf  /( h2aUa ) .  If  the  calculated  value  of  hcWVf  l(h2aUa )  is  larger  than  1,  then  p  =  1,  and  all 
material  will  fall  into  the  channel. 

If  /^-percent  of  the  suspended  material  at  a  given  instant  is  trapped  in  the  channel,  then  (1  -p) 
percent  bypasses  the  channel  by  being  advected  over  it.  Kraus  and  Larson  (2001)  describes  how 
material  may  also  be  resuspended  and  transported  out  of  a  channel,  thereby  increasing  the 
amount  bypassed. 

The  depth-averaged  current  velocity  perpendicular  to  the  channel  should  be  employed  for 
predicting  the  trapping  ratio.  If  the  current  approaches  the  channel  at  an  angle  0  with  respect  to 
the  orientation  of  the  channel  (i.e.,  for  a  perpendicular  current  0  =  90  deg),  the  current  should  be 
divided  by  sin  0  before  it  inserted  as  Ua  in  Equation  7.  This  result  indicates  that  the  trapping 
ratio  increases  because  the  distance  traveled  for  the  sediment  particles  across  the  channel 
becomes  longer  (or  equivalently,  the  component  of  the  incident  current  transporting  the  particles 
across  the  channel  becomes  smaller).  If  there  is  a  unifonn  current  along  the  channel  axis  during 
its  crossing  of  the  channel,  the  end  result  of  consideration  of  the  processes  is  that  there  is  no 
change  in  the  trapping  ratio.  Although,  the  distance  the  sediment  particles  travel  increases,  so 
does  the  speed,  and  the  two  processes  cancel. 


Figure  3.  Percentage  of  sediment  falling  into  the  channel 


EXAMPLES:  The  preceding  development  provides  a  means  of  estimating  the  percentage  of 
suspended  material  that  may  fall  into  a  channel  (or  which  may  bypass  a  channel)  based  on 
sediment  fall  speed,  geometry  of  the  channel,  and  the  speed  of  the  upstream  flow  crossing  the 
channel.  Here,  two  examples  are  given  to  illustrate  applications  of  this  information. 
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Example  1 :  Material  trapped  by  or  crossing  a  channel  during  (a)  typical  surf  conditions, 
(b)  and  during  a  storm. 

Given:  Channel  width  W  =  30  m;  channel  depth  hc  =  4.5  m;  average  ambient  depth  ha  =  3  m; 
sediment  grain  size  =  0.2  mm.  During  typical  surf  conditions,  Ua  =0.1  m/s  and  nearshore  wave 
height  is  1  m;  during  stronger  storms,  Ua  =1  m/s  and  nearshore  wave  height  is  4  m. 

Find:  The  relative  percentages  of  suspended  material  deposited  into  and  bypassing  the  channel 
for  that  portion  of  the  channel  located  in  the  surf  zone. 

Solution:  (a)  Typical  wave  condition:  Waves  of  1-m  height  will  break  in  water  depth  of  about 
1.3  m,  or  in  depths  well  shoreward  of  the  area  of  interest  at  3-m  depth  and  seaward  of  the  jetties. 
However,  because  of  the  wave  current  interaction  and  the  presence  of  the  ebb  shoal  in  the 
vicinity  of  the  channel,  wave  breaking  is  observed  to  occur  intermittently,  suspending  sediment. 
Therefore,  we  proceed.  Consulting  Table  1,  for  0.2-mm  sand  at  20-deg  C,  we  have  V/  = 
0.025  m/s,  giving: 

Fr  =  Vf  /  Jgha  =  0.025  /  V9.8x3  =  0.0046 , 


and 


hc WVf  /(h2aUa)  =  4.5  x  30  x  0.024 /(9 x  0. 1)  =  3 .75 

Calculating  Equation  7  or  entering  Figure  2,  one  finds  that  p  >1  for  any  value  of  Fr,  so  p  = 
100  percent,  and  all  the  material  falls  into  the  channel  for  this  condition. 

(b)  Stonn  condition:  Waves  of  4-m  height  will  be  breaking  all  around  the  channel.  For  this 
situation: 

Fr  remains  as  previously  described,  0.0046,  and 
hcm 7  /(h2aUa)  =  4.5  x  30  x  0.025  /( 9  x  1)  =  0.375 

Calculating  Equation  7  or  entering  Figure  2  for  the  curve  Fr  =  0.005,  p  =  46  percent,  which 
means  that  1  -p  =  54  percent  of  the  material  passes  the  channel.  It  is  again  noted  that  not  all 
suspended  material  that  falls  into  the  channel  will  remain  there.  Some  portion  will  be 
resuspended  and  carried  out  of  the  channel,  thereby  increasing  the  percent  bypassed. 

Example  2:  Determine  the  possibility  that  a  sediment  plume  viewed  in  an  aerial 
photograph  consists  of  sand  bypassing  an  inlet. 

Given:  An  aerial  photograph  at  an  inlet  indicates  that  a  sediment  cloud  or  plume,  evident  by 
brown  color  observed  on  the  background  of  green  and  white  turbulent  water,  sometimes  appears 
in  the  surf  zone  under  typical  wave  conditions  (wave  height  less  than  1  m)  that  follow  after  a 
storm.  The  plume  bypasses  both  jetties  of  the  inlet  entrance  and  continues  alongshore  for  a 
distance  of  200  m  past  the  updrift  jetty,  where  the  plume  originates.  The  water  depth  near  the 
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ends  of  the  jetties  is  about  4  m,  and  the  channel  depth  is  6  m.  The  grain  size  at  this  beach  is 
relatively  fine,  0.17  min,  and  has  a  fall  speed  of  0.02  m/s.  The  typical  longshore  current  speed 
under  fair-weather  wave  conditions  is  estimated  to  be  0.2  m/s,  a  fairly  strong  current  because  of 
the  strong  longshore  component  of  the  wind  at  this  site. 

Find:  If  it  is  assumed  that  the  sediment  plume  can  be  observed  half  way  through  the  water 
column,  is  it  possible  that  the  downdrift  portion  of  the  plume  consists  primarily  of  0.17-mm 
sand? 

Solution:  Under  typical  wave  conditions,  waves  cannot  break  in  the  channel,  and  resuspension 
will  be  minor.  Therefore,  we  assume  the  suspended  sand  originates  from  the  surf  zone  and  is 
carried  by  the  rip  current  along  the  updrift  jetty  to  be  captured  again  by  the  wind-generated 
longshore  current  and  transported  past  the  jetties.  We  take  Az  =2  m  to  be  the  maximum  elevation 
in  the  water  column  that  can  be  observed  downdrift  and  assume  some  0.17-mm  particles  are 
present  at  the  water  surface  at  the  updrift  (starting)  end  of  the  plume.  Then  by  Equation  6, 
interpreting  the  distance  W  as  the  maximum  length  of  travel, 

W  =  —  —  Az  =  --^-2  =  20  m 
hc  Vf  4  0.02 

In  this  calculation,  we  set  the  ratio  hjhc  =  1  to  maximize  the  distance.  The  fine  sand  suspended 
to  the  water  surface  and  carried  by  a  0.2-m/s  current  can  only  travel  20  m  before  dropping  below 
2  m  in  the  water  column.  Therefore,  it  is  concluded  that  the  surface  portion  of  the  plume 
observed  in  the  photograph  must  consist  of  fine  material  such  as  clay  that  has  small  fall  speed 
and  of  organic  material  that  floats.  Because  such  plumes  can  appear  after  a  storm,  it  may  be  that 
the  material  consists  of  clay  and  organic  particles  washed  into  the  ocean  during  the  runoff  caused 
by  the  storm. 
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Center  (ERDC),  Coastal  and  Hydraulics  Laboratory.  The  research  was  jointly  supported  by  the 
Coastal  Inlets  Research  Program,  Inlet  Channels  and  Adjacent  Shorelines  Work  Unit,  and  by  the 
Coastal  Navigation  and  Sedimentation  Program,  Diagnostic  Modeling  System  Work  Unit. 
Questions  about  this  CHETN  can  be  addressed  to  Dr.  Nicholas  C.  Kraus  (601-634-2016, 
Fax  601-634-3080,  e-mail:  KrausN@wes.army.mil). 

This  CHETN  should  be  cited  as  follows: 

Kraus,  N.  C.,  and  Larson,  M.  (2001).  “Estimation  of  suspended  sediment  trapping 
ratio  for  channel  infilling  and  bypassing,”  ERDC/CHL  CHETN-IV-34,  U.S.  Army 
Engineer  Research  and  Development  Center,  Vicksburg,  MS. 

(, http://chl.wes.army.mil/library/publications/chetn ) 
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